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Abstract 

> 

^ \ We study the possibility that cosmological density perturbations are generated by the inhomogeneous 
decay of right-handed neutrinos. This will occur if a scalar field whose fluctuations are created during 
inflation is coupled to the neutrino sector. Robust predictions of the model are a detectable level 
of non-Gaussianity and, if standard leptogenesis is the source of the baryon asymmetry, a baryon 
isocurvature perturbations at the level of the present experimental constraints. 

1 Introduction 

Recently a lot of attention has been devoted to models of the early Universe where the leading 
source of density perturbations is not the slowly rolling inflaton but an additional light scalar. The 
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quantum fluctuations of this additional field generated during inflation can be converted in different 
ways into adiabatic density perturbations. In the curvaton scenario [T] the scalar starts oscillating 
when the Hubble scale drops below its mass and its energy density becomes quite relevant before 
decay. As in different parts of the Universe the initial conditions for the oscillation are slightly 
different this gives rise to density perturbations. Another possibility is that the scalar fluctuations 
modulate some parameter relevant for the cosmological evolution at some stage. For example it can 
modulate the reheating efficiency of the inflaton or the mass and interactions of a particle which 
becomes dominant before decaying [21IH1II]- 

These models are theoretically interesting: given the ubiquity of scalar fields in extensions of the 
Standard Model, it is easy to imagine that the inflaton is not the only relevant degree of freedom 
in the early phases of the Universe. In particular, in String Theory dimensionless couplings are 
functions of the scalars describing the compactification. Moreover this additional source of density 
perturbations can allow to circumvent the necessity of slow-roll altogether |Sj . If one takes the point 
of view that the inflaton parameters are in some way scanned over a landscape of models, extreme 
values for the amplitude of density perturbations (either much smaller or much bigger than the 
observed 10 -5 ) seem preferred (El El- This is easy to understand as changing the inflaton parameters 



1 



one also changes the number of e-folds of inflation, so that the final volume exponentially depends 
on the parameters. If taken seriously, this can be a further motivation for looking at alternative 
sources of perturbations. 

However the main reason behind the study of these models is that they give experimental sig- 
natures which can allow to distinguish this source of perturbations from the standard inflaton fluc- 
tuations. There are two basic observables which are relevant. First, in general these models predict 
a higher level of non-Gaussianity with respect to the standard slow-roll scenario, with a well de- 
fined momentum dependence. Second, given the presence of a second field, they can give rise to an 
isocurvature contribution. 

In an attempt to connect these ideas with realistic particle physics scenarios, in this paper we 
will study in detail a well defined and simple model. We will assume that the additional light scalar 
is coupled to right-handed (RH) neutrinos, so that their mass and decay rate are slightly different 
in different regions of the Universe. This will give rise to density perturbations. As we will see 
the production of density perturbations requires departure from thermal equilibrium. Therefore 
RH neutrinos appear as the most simple and minimal candidate: being SM gauge singlets, they 
decay out of equilibrium for sufficiently small Yukawa couplings. Fluctuations in the RH neutrino 
parameters will also modulate the leptogenesis process (which we assume to be the source of baryon 
asymmetry), thus giving rise to a baryon isocurvature contribution, correlated with the adiabatic 
one. We will see that the level of non-Gaussianity and the isocurvature contribution put severe 
constraints on the scenario. Although model dependent, the amplitude of the baryon isocurvature 
mode is generically predicted to be somewhat bigger than the present limits. Given this tension, 
the detection of baryon isocurvature should be around the corner: the model will be ruled out if the 
next generation experiments do not detect anything. 

There are various reasons why baryon isocurvature is interesting, especially in this class of models. 
First of all, given that baryon number is conserved 1 , it is much easier to produce baryon isocurvature 
than for example dark matter or neutrino isocurvature. These last two are in fact erased by thermal 
equilibrium. 

A second point is that both the generation of the baryon asymmetry and the production of density 
perturbations through the modulation of some parameter are processes which require departure 
from thermal equilibrium. This is well known for baryogenesis, let us see how it works for the other 
process. Density perturbations are produced in these models because the cosmological evolution is 
different in distant regions of the Universe as some relevant parameter is made space dependent by 
the fluctuations of the light scalar. In thermal equilibrium entropy in a comoving box is constant. 
This gives a unique relationship between the scale factor and the temperature of the plasma and 
thus it forbids the creation of density perturbations. All the explicit models are in fact based 
on out of equilibrium processes: reheating after inflation, inhomogeneous dominance of a massive 
particle (which will be the example of interest here) and one could imagine further examples like the 
modulation of a phase transition. Now if this out of equilibrium process is also the source of baryon 
asymmetry then one expects the generation of baryon isocurvature at some level. For example it 
is difficult to imagine that the reheating process is completely baryon symmetric 2 0, so that some 

1 Given the presence of sphalerons at high temperature what we really mean is that B — L is conserved. 
2 Notice that as the thermal phase of the Universe begins after the decay of the inflaton, we are not forced 
to look for additional sources of the baryon asymmetry. The simplest scenario is that the thermal phase 
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baryon isocurvature will be generated in a inhomogeneous reheating scenario. The possibility of a 
baryon isocurvature is well known also in the context of curvaton models 

Finally the baryon isocurvature experimental limits are entering now in a very interesting range 
from the point of view of the models we are studying. From what we said, we expect the baryon 
isocurvature contribution 5{ub / s) / (ns / s) (ng is the baryon density and s the entropy density) 
to be comparable to the adiabatic component which can be parametrized by the usual variable Q: 
S(ns I ' s) I '{jib 1 1 s) ~ (. The most recent global data analysis jTUJ (see also puts a constraint 3 

< 0.7 at 2cr . (1) 

Thus we are really in a very interesting range! 

In the following Section we will calculate the perturbations created by the modulation of the 
RH neutrino decay, focusing on two possible signatures of the model: non-Gaussianity and baryon 
isocurvature. We start assuming that the fluctuations in the underlying scalar field are frozen, so 
that their only effect is to change the neutrino parameters. We will study the scalar evolution in 
Section El before drawing our conclusions in Section 0J 



S{ub/s) 



n B /s 



2 Generation of density perturbations 

In this section we want to study the production of density perturbations through the inhomogeneous 
dominance and decay of right-handed neutrinos. Fluctuations in the RH neutrino parameters comes 
from the coupling with a scalar eft, whose super-horizon perturbations are generated during inflation 4 . 
We stick to the simplest framework of Standard Model plus three hierarchical RH neutrinos. We 
concentrate for now on the decay of the lightest one and on adiabatic perturbations, postponing the 
discussion about fluctuations in the ratio ub/s. We assume that this mechanism is the leading source 
of density perturbations; we will see that this implies that the RH neutrino must go rather strongly 
out of equilibrium. If this is not the case its density remains a small fraction of the plasma which 
contains ~ 100 relativistic species; therefore to match the observed ~ 10~ 5 level of perturbations 
much bigger fluctuations in the decay width 5T /T and mass 5M/M are required. This would imply, 
as we will see explicitly below, that second order corrections are quite big and incompatible with 
the present limits on non-Gaussianities. 

The ratio T/H controls whether or not the RH neutrino decays are in equilibrium. In particular 
it is useful to study the quantity 

T(T = 0) (AtA)n-M / 1/2 M 2 2 tt 3 / 2 x 
H(T = M) 8tt 1 V* M P ^45 / ' 

where Mp = G^^ 2 ~ 1.22 x 10 19 GeV and g* is the number of relativistic degrees of freedom in 
the plasma = 106.75 for the Standard Model). The mass of the lightest RH neutrino is M and 

just starts from an asymmetric initial condition due to a baryon asymmetric reheating. Obviously it is still 
interesting to try to relate the observed baryon abundance to known physics, as in the case of leptogenesis. 

3 The constraint holds assuming that the perturbations are completely correlated or anti-correlated and 
have the same scale dependence. This applies to our model as both kinds of perturbations come from the 
same scalar </>. 

4 Additional interactions between cf> and the Standard Model are possible but irrelevant in our discussion. 
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the matrix A gives the Yukawa couplings of neutrinos with the Higgs. The ratio can be rewritten 
introducing the Higgs VEV at low energy v {v ~ 174 GeV) in the form 

r(T = 0) mx „ _(AU)n^ 

— where mi = — . (6) 



H{T = M) m* ' M 

The quantity mi (the usual decay parameter used in leptogenesis) roughly corresponds to the con- 
tribution of the lightest RH neutrino to the left-handed neutrino mass matrix. On the other hand 
m* is a fixed energy scale given by 

For mi <C m* RH neutrinos are out of equilibrium and decay at a temperature much smaller than 
their mass. If the decay happens very late neutrinos dominate the Universe before decaying. This 
will happen after the Universe has expanded by a factor since T ~ M. As in radiation dominance 
H oc a~ 2 we need rh\ <C rn*/g 2 . 

We are interested in the density perturbations induced by the inhomogeneous dominance and 
decay of RH neutrinos; we assume that RH neutrinos start from a thermal distribution for T 3> M. 
In the absence of additional sources of perturbations, at temperature T 3> M the Universe is 
unperturbed and the surfaces of constant temperature are flat. Different regions of the Universe 
will be characterized by different values of the scalar eft, so that the mass and decay width of the 
RH neutrino will be different. As a consequence the neutrino will become non-relativistic and 
decay at slightly different temperatures in different parts of the Universe. If we take two reference 
temperatures, one above the RH neutrino mass (let us call it T^h) and one much below, after the 
neutrino decay (Ti ow ), the ratio of the scale factors between these two temperatures a(Ti ow )/a(Thi g h) 
will not be exactly the same in different parts of the Universe. Thus if the surface of constant Thigh 
is flat, the one of Ti ow will be curved. Modes of cosmological interest are huge compared to the 
horizon at the temperature we are looking at. In this limit adiabatic scalar perturbations can be 
expressed, in coordinates in which the surfaces of constant time are of constant temperature, as 

ds 2 = -dt 2 + e 2 ^a(t) 2 dx 2 . (5) 

We have neglected terms subleading in the expansion k/(aH), where k is the comoving wavevector 
of the perturbation. In this limit the variable £ is constant in time at any order of perturbation 
theory. The reason for this is quite clear: £ just describes a locally unobservable rescaling of the 
spatial coordinates ^21 QH| • From the discussion above we conclude that perturbations induced by 
the RH neutrino can be obtained by the equality 

e^) = ^4(M,r). (6) 

Spatial fluctuations of M and T induce fluctuations in £. In this way the problem of calculating den- 
sity perturbations is reduced to evaluate the right hand side of the equation above in an unperturbed 
Universe as a function of the parameters describing the RH neutrino. 

Let us start from the simple case in which the RH neutrino dominates before decaying. As we 
discussed this will happen for rh\ <C rh* jg 2 . In this case we can divide the evolution of the Universe 
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around this epoch into three periods with different equation of state. (1) At temperatures T 3> M/g* 
the Universe is radiation dominated. (2) The Universe is dominated by the RH neutrino (matter 
dominance) from T ~ M/g* until its decay at H ~ T. (3) After the neutrino decays we are back in 
radiation dominance. It is easy to calculate the expansion in the three stages 

a(T ~ M/g*) 



a (^high) 



oc M (7) 



a(H~T) ( H(T ~ M/g* 



a(T ~ M/g* 



OC 



2/3 



±-^n±L \ ocM 4 / 3 r- 2 / 3 (8) 



a(T low ) / T 



OC 



1/2 



a(ff ~ r) 

We can put these results together to get j3j 

a(Ilow) ,,l/3p-l/6 _ ~-l/6 



oc 



M l/3 p-1/6 a ~- 1/6 for ~ i < ~*y^2 _ ^ 



« (Thigh) 

Note that we are not interested in the proportionality factor as it corresponds to a space indepen- 
dent rescaling of a. In the limit we studied the result just depends on rh\ and not on M and T 
independently. It is easy to see that this holds in general. If we change M keeping fh\ fixed we are 
moving the window of temperature where the equation of state deviates significantly from radiation 
dominance, but this does not change the expansion between the reference temperatures Thigh an d 
Tj ow . In fact, once we fix fh\, the scale factor becomes a function of the ratio T/M only and not of 
T and M independently; as both at high and at low temperature we are in radiation dominance the 
dependence on M simplifies in the ratio a(T\ ow / 'M) / 'a(Thi g h / 'M) . Thus relation Q simplifies to 

OUhighJ 

Let us now relax the inequality rh\ <C fh*/g#- For larger values of rh\ the RH neutrino does not 
dominate the plasma before decay; increasing rh\ the deviation from simple radiation dominance 
becomes smaller and smaller until eventually the ratio of scale factors a(?] ow )/a(Thigh) approaches 
the constant value Thi g h/T| ow . The solution for a generic rh\ can be easily obtained numerically by 
solving the coupled equations describing the decay of the RH neutrino fluid into radiation 

Pl + 4H Pj = T PN (12) 
p N + 3(l + w N (T/M))Hp N = -T PN (13) 

H 2 = ^liPN + P,)- (14) 

We are assuming that the RH neutrino has a thermal population at very high temperature so that 
its energy density pn starts as 1/g* of the one in the plasma p 7 . The pressure to energy density 
ratio wn = Pn / Pn of the neutrinos drops from 1/3 to as they become non-relativistic at T ~ M. 
We will see that the decay must occur sufficiently out of equilibrium to be consistent with the limits 
on non-Gaussianities, so it is self-consistent to neglect in the equations the inverse decay of the RH 
neutrino. We also neglect AL = 1 reactions which give contributions suppressed by a loop factor 
and thermal corrections, which are small as the decay occurs at T <C M |14j . In Fig. ^we plot the 
ratio a(?] ow )/a(Thigh) as a function of the parameter rh\. For rh\ <C fn*/g 2 we recover the analytic 
result of eq. (|lfl|) : while for much bigger rh\ we approach the radiation dominance result. 
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Log(a(Ti 0W )/a(T high ) Ti ow /T high ) 




Figure 1: Expansion between two fixed temperatures as a function of mi. 

2.1 Non- Gauss ianity 

Now that we have calculated the expansion ratio as a function of fhi, we are ready to evaluate 
the amount of non-Gaussianity. Relation Ijlljl between rh\ and £ is non-linear and this is the 
source of non-Gaussianity. Let us assume for the moment that the ratio Srhi/rhi is distributed as a 
Gaussian. Expression gives £ in terms of rh\ at the same point in space. In this case, the form 
of non-Gaussianity induced by second order terms is usually called local 5 and parametrized by the 
parameter /nl defined by 

C(x) = Q(x) - ^h L (C g (x) 2 - (C 2 g )) , (15) 

where Q g is a Gaussian random variable. 

From the function £ = f (log fhi/fh*) shown in Fig. ^ we have at second order 

C(S) = /'^(: ? ) + i(r-/')(ai(i>)) 2 , (16) 
mi 2 \ mi J 

where we have subtracted an irrelevant constant value in Q. From the linear term we get the 
relationship between the spectrum of Sfhi/rhi and the one of £, while the second order term gives 

/nl = ~ - ■ (17) 

As we expected this expression becomes very large for large rh\. The RH neutrino becomes less and 
less relevant at decay, so that keeping fixed the observed value 10~ 5 for density perturbations, we 

5 The shape dependence of the 3-point function in different models has been studied in [T^lUT)] . 
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have to correspondingly increase the fluctuations 5fh\/fh\. This implies that second order correc- 
tions become more and more relevant. This is an example of a quite generic correlation between 
the "inefficiency" of the mechanism producing density perturbations and the level of non-Gaussianity 
|17j . The non-linearity parameter /nl will be of the order of the inverse of the fraction of energy 
density in the RH neutrino at decay: 

|/nl|~3*\/I|. (18) 
V m* 

In Fig. [2 we plot /nl as a function of rh\ from eq. (|17[). 




In the discussion we assumed that the ratio 8fh\/fhi is distributed as a Gaussian. This will not 
be the case in general; the dependence rh\{<f>) will induce additional non-linearities similarly to what 
happened in the relation ((mi). It is easy to realize that also these non-linearities will be enhanced 
by the "inefficiency" of the mechanism when the RH neutrino does not dominate. The coefficient 
/' in eq. becomes small in this limit and second order terms in the expansion of rhi(4>) give a 
contribution to /nl enhanced by / . We conclude that although the plot of Fig. El should not be 
taken as a sharp prediction we do not expect, barring unlikely cancellations, a smaller value for /nl 
( 6 )- 

In addition to the non-linearities we discussed the statistics of (f> itself may be non-Gaussian as a 
consequence of its self-interactions. The leading effect will be given by the non-linear superhorizon 
evolution, which gives a non-Gaussianity of the form 1)151) [17] . Also this contribution is quite model 
dependent and again enhanced by the "inefficiency" / . 

The most stringent experimental limits on the parameter /nl come from the analysis of WMAP 
data [El EH- The allowed 2a range —27 < /nl < 121 can be converted into an approximate (given 

6 Note that even the negative sign for /nl in Fig. [21 can be flipped by second order effects from fh\ (</>). 
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the unknown contributions discussed above) bound on mi 

mi < 1(T 3 • m* ~ HT 6 eV . (19) 

This means that the RH neutrino Yukawa interactions are strongly out of equilibrium at decay. This 
implies that the same interactions cannot bring the neutrinos in thermal equilibrium at T > M; we 
will come back to this point in the following Section. 

Notice that in general inequality Q19JI has no direct implication for the low energy neutrino 
parameters as the neutrino parameters relevant for the production of density perturbations are 
different from the ones measured at low energy. Indeed the scalar starts displaced from the origin 
where it will be eventually stabilized. If one assumes that the variation of parameters is not very 
large, the inequality above implies that one of eigenvalues of the left-handed neutrino mass matrix 
is smaller than 10 -6 eV. In fact it is easy to prove that m± cannot be smaller than all the three 
neutrino mass eigenvalues |2U| . From the constraints we have from neutrino oscillation experiments, 
this state would be the lightest of the three, with a resulting hierarchical spectrum. 



2.2 Baryon isocurvature 

We now come to the discussion of the baryon isocurvature contribution. We will see that this puts 
additional severe bounds on the model. The ratio between baryon and entropy density produced by 
RH neutrino decay after reprocessing by sphaleron transitions can be written as ^1] 

J = -f 9 e NlV (fn 1 )^f(T»M). (20) 



The factor ejvi is the CP-asymmetry parameter in Ni decay which is given by 

1 hn(AtA)j t Jm%. 
Svr^ (AtA)n {m 2 Ni 

(1 + x) log 



^ ~ ■ (22) 
2yjx 



x 

The ratio n^/s is calculated in thermal equilibrium at high temperature and it is given by n^/s = 
135C(3)/(4vr 4 ^). Finally 7] is an efficiency parameter 7 which takes into account the wash-out of 
the produced baryon asymmetry (which is relevant if the RH neutrino decays close to equilibrium 
rh\ > fh*) and the entropy released in the decay (which is relevant when the RH neutrino becomes 
dominant in the Universe fh\ < m* /gl). 

Let us start assuming that the CP violating parameter ejvi does not depend on the scalar <p so 
that it is constant in space. Given the constraints we discussed from non-Gaussianities, the decay of 
the RH neutrino must be strongly out of equilibrium so that the wash-out of the produced baryon 
asymmetry can be neglected. Every neutrino becomes (on average) a constant fraction —28/79 • 
of a baryon, independently of rh\. This implies that the surfaces of constant baryon density remain 



7 The efficiency parameter depends solely on mi, unless the mass of the RH neutrino is very large: M > 
10 13 GeV PI|. 
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flat while, as we discussed, the surfaces of constant temperature become curved after the RH neutrino 
decay. This mismatch gives a baryon isocurvature component 

S(n B /s) = _5s = _ 3 6T = _ 3C 
tib/s s T 

This isocurvature component is too large and it is ruled out by the present limits in Additional 
contributions will come from fluctuations in the CP violating parameter ejvi ■ if the scalar cj) enters in 
the mass term and/or in the Yukawa interactions of the RH neutrinos this will induce fluctuations 
in ejvi • Thus we have 

6(n B /s) _ 3 | 5e Nl /e Nl 
n B /s ( 

Roughly we expect also this second contribution to give an order one isocurvature component. From 
the explicit expression for ejVi eq. (|2"Tj) . we see that the precise number is very model dependent 
as it depends on how the scalar <p enters in the different entries of the Yukawa matrix and in the 
right-handed neutrino masses. 

We conclude that a moderate cancellation between the two terms above is required to be consis- 
tent with the present limits. From a more optimistic point of view we can say that detection of an 
isocurvature component should be around the corner, as it is naturally expected to be of order one 
in these models. 

What happens if we consider not only the lightest RH neutrino but also the other two (assuming 
that they all start with thermal abundance at high temperature)? We can study two interesting 
cases. In the case all the three neutrinos are not in thermal equilibrium when they become non- 
relativistic (m; <C rh* for all i), then the situation is pretty similar to the case studied below with 
just a single RH neutrino. All the three neutrinos will contribute to the generation of £ (we can 
take the reference temperatures Thigh and Ti ow to encompass all the neutrinos) and to the baryon 
asymmetry. Note that as all the neutrinos are out of equilibrium the baryon asymmetry generated 
by a RH neutrino is not washed out by the lighter ones. Again if we take all the CP violating 
parameters ej as constant we get an isocurvature component 3 times bigger than the adiabatic one 
as in eq. (|23[), Thus again a moderate cancellation between this contribution and the one generated 
by Sti/ei is required. 

A quite different situation is the case in which the lightest RH neutrino N% is close to equilibrium 
when it becomes non-relativistic. As we discussed its contribution to £ will be small in this case 
and can be neglected. As N% is in equilibrium after the heaviest neutrinos decay, it will wash- 
out the lepton asymmetry already produced. Although this wash-out effect depends on the flavor 
structure and can be not entirely efficient |2H I22| , we can assume for simplicity that all the baryon 
asymmetry is generated by the lightest RH neutrino. Thus now the two mechanisms are separated: 
fluctuations in the parameters of the heavy generate £ while fluctuations in the decay of N% 
give baryon isocurvature. In the limit in which the scalar <j> is decoupled from N% {i.e. it does not 
change neither its mass nor its couplings) we have no isocurvature. As N\ is close to equilibrium 
the wash-out of the produced lepton asymmetry cannot be neglected. A good parameterization of 
the efficiency parameter fy(mi) in the regime rh\ S> rh* is given by r] oc mj -1 ' 16 |14j . Fluctuations 
in this efficiency parameter and in will give an isocurvature component. Thus, although in 
this case we do not have to rely on a cancellation to be compatible with data, we still expect the 
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isocurvature contribution to be of order one, unless for some particular reason the coupling of (j) with 
N\ is suppressed. 

More complicated intermediate scenarios are possible but, without further specifying the coupling 
of (/>, the generic conclusion we can draw is that the isocurvature contribution should be around the 
(experimental) corner. 

Up to now we assumed that, at high temperatures, RH neutrinos are in thermal equilibrium. On 
the other hand we need that, at least for one of the neutrinos, the decay parameter fh is sufficiently 
small to allow the RH neutrino to become a significant fraction of the plasma before decay. This also 
implies that Yukawa interactions are not fast enough to bring the neutrino in thermal equilibrium 
starting from a generic initial condition. Additional interactions are required to produce neutrinos 
at high temperature; for example gauge interactions mediated by a heavy GUT gauge boson 23 . 
If there are no additional interactions we expect RH neutrinos to get produced at some level by the 
reheating process, either perturbatively or by their coupling with the coherent infiaton oscillations 
|24j . In this case the initial RH neutrino density is model dependent and we do not even know their 
pressure to energy ratio [25]. The analysis above obviously does not apply to this case, although we 
still expect a significant baryon isocurvature component as the generation of density perturbations 
and baryogenesis are unified in the same process. We can think about this non-thermal production of 
RH neutrinos as part of the reheating process, before getting to a completely thermal state. From this 
point of view the modulated decay of RH neutrinos becomes a modulated reheating scenario [2111]. 
This is a nice example of a more general point: if the baryon asymmetry is produced at reheating 
(which is the simplest mechanism for baryogenesis), we expect that models with modulated reheating 
will give a baryon component. 



So far we concentrated on the effects of the modulation of the RH neutrino parameters induced by 
the fluctuations of a scalar <p. Another aspect one should consider is the effect of the scalar couplings 
with the neutrino sector on the dynamics of the scalar itself. We have tacitly assumed that the scalar 
does not evolve in time around the epoch when RH neutrinos become non-relativistic, so that its 
only effect is to make the parameters space dependent. We can now ask if this is a good assumption. 
Let us start assuming that (j) enters only in the RH neutrino mass matrix with a coupling (neglecting 
all the flavor structure) 



Without loss of generality we choose to normalize the scalar with Mp in the function M. The scalar 
is coupled to the plasma of RH neutrinos; assuming that they start with a thermal distribution we 
have (NN) ~ T 3 , so that the evolution of <fi will satisfy 



3 The evolution of the scalar <j) 




(25) 



<p + m<p + 



M' 



T 3 = . 



(26) 



M P 
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We are assuming that the bare potential for 0, V(cf>), is sufficiently flat to be irrelevant at the epoch 
we are looking at 8 . The variation of <f> during an Hubble time will thus be of order 

A , M' MT 3 , . 

AS — T Mp . (27) 

As we need the RH neutrino to become dominant to be a sufficiently efficient source of density 
perturbations, during dominance we have MT 3 / \H 2 Mp) ~ 1 which gives 

Acpy^Mp. (28) 
On the other hand the induced level of density perturbations will be at most of order 

M> 6(/) M' //inflation 

s M M P M M P y J 

From the fact that we did not observe gravitational waves created during inflation, we have an 
upper bound on inflation which, given the observed value of £, gives a lower bound on M'/M of 
order M'/M > 10 ( 9 ). This implies that the scalar moves by more than Mp during the epoch we 
are interested in! Even taking the natural range for <j> of order Mp, this motion will induce a very 
big variation in the neutrino parameters. This motion makes all the analysis more complicated and 
dependent on the stabilization potential of the scalar. 

The situation is much more under control in the case the scalar enters only in the Yukawa 
couplings of RH neutrinos 

X(dj/M P )LHN . (30) 
In this case, as (LHN) ~ AT 4 , we have an equation of motion 

4> + 3H<}> + ^-T 4 = . (31) 
The variation of dj in one Hubble time is thus 

A' A 2 T 4 

a *~JhH4 Mf - <32) 

while the expression for density perturbations is 

\/ IT 

/\ -TC inflation /oo\ 

C ~A^p-' (33) 

We see that, contrary to the previous case of mass variation, for small enough A the displacement 
of the scalar is negligible without suppression of the produced density perturbations. We conclude 
that this scenario with fluctuating Yukawa couplings is favored with respect to mass fluctuations. 



8 This requirement limits the amount of non-Gaussianity coming from the self- interactions of <f> |17| . 

9 Notice that we are not using the constraint 1201 that the level of density perturbations created by 
our mechanism is dominant with respect to the standard contribution from the inflaton Q ~ H / \Mp\fe). 
This constraint (which would give a lower bound on M'/M a bit stronger than the one above) could be 
circumvented in models of inflation without slow-roll jS] ■ On the other hand the gravitational wave constraint 
discussed in the text is completely general. 
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We have discussed the evolution of (ft at the time density perturbations are generated, but what 
happens to it after that? Eventually the mass term in the scalar potential V((ft) will get larger than 
the Hubble rate H. At this point the scalar starts oscillating around the minimum and the energy 
stored in these oscillations red-shifts as a -3 , i.e. slower than radiation. If the displacement of the 
scalar from the minimum is of order Mp, the oscillations start dominating soon after > H. 
In this case the scalar would act as a curvaton PQ: everything created before gets redshifted away 
(including the generated baryon asymmetry) and, as the oscillations have a slightly different initial 
conditions in different regions, we have an additional source of density perturbations, which adds to 
the inhomogeneous RH neutrino decay. 

Whether the Universe becomes dominated by the scalar oscillations or not is, however, very 
model dependent. First, its amplitude of oscillation could be much smaller than Mp either because 
of initial conditions or because the typical scale of variation for (ft is much smaller than Mp. This 
happens for instance for an axion-like scalar, where the typical range of variation is the decay 
parameter /, in general much smaller than Mp. Second, its decay into some hidden sector or into 
the visible one can be fast enough to avoid the domination. In the second case one must be careful 
that the additional interactions are not so strong to induce a huge displacement of (ft before RH 
neutrino decay. Given our complete freedom (ignorance) about the coupling of (ft, it is easy to check 
that this can be done. 

Another interesting possibility about the fate of (ft would be that the scalar is still light in the 
present Universe and possibly observable |27j . 

4 Conclusions 

In this paper we studied a concrete realization of the inhomogeneous decay scenario (3] which is 
based on standard thermal leptogenesis. We considered the simplest possibility: the SM extended 
with 3 hierarchical RH neutrinos. Both curvature perturbations and the baryon asymmetry are 
produced through the out-of-equilibrium decay of RH neutrinos whose parameters (mass and Yukawa 
couplings) are controlled by a light scalar (ft. 

The most interesting feature of the model is that a baryon isocurvature mode is produced which 
is correlated with the adiabatic one. Though the amplitude of this mode is model-dependent, we 
generically expect that detection should be very close. Focusing on the lightest RH neutrino, the 
dynamics is completely described by a single parameter rh\. It sets the departure from equilibrium 
at decay, which controls the efficiency of both leptogenesis and the production of perturbations. 
The present limits on non-Gaussianity constrain rh\ to be in the weak wash-out (strongly out-of- 
equilibrium) regime. 

We also studied the evolution of the scalar field responsible for perturbations. The couplings of 
the scalar with the neutrino sector induce a motion of the scalar field itself, which in turn induces 
variations in the neutrino parameters. This puts additional constraints on the model. In the case 
where (ft enters only in the mass of the RH neutrinos, (ft is pulled over many Planck masses during 
RH neutrino domination. This possibility is therefore disfavored. On the other hand, if (ft enters 
only in the Yukawa couplings A, the pull is suppressed by A 2 so that it can be made sufficiently 
small. 

Many generalizations of the model are possible, the simplest being an embedding in supersymme- 
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try. In this case the right-handed sneutrino is an additional source of baryon asymmetry. Moreover 
the possibility of using a sneutrino direction as inflaton or curvaton has been thoroughly studied 
[SI El- This opens a broad range of model building possibilities. 
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